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.  __  ABSTRACT 

This  thesis  is  a  continuation  of  a  series  of  studies  on  automatic 
throttle  compensation  systems  begun  in  1963  by  Lieutenant  Commander 
G.  R,   Bell.  £l]   Previous  investigations  have  been  concerned  with  flight 
in  a  wings  level  attitude.     This  study  will  consider  the  effects  of  turning 
flight  on  the  system.     The  system  under  consideration  here  is  the  Special- 
ties,  Incorporated,  Automatic  Power  Compensator,    (APC),  an  automatic 
throttle  in  use  in  the  F-8  and  F-4  aircraft.     It  utilizes  feedbacks  of  angle 
of  attack  and  normal  acceleration.     It  was  found  important  to  differentiate 
between  spurious  perturbations  and  those  perturbations  from  steady  state 
introduced  due  to  the  aircraft  entering  a  turn. 

This  study  was  conducted  by  Lieutenant  Richard  D.   Hartman  at  the 
United  States  Naval  Postgraduate  School,   Monterey,  California.     The 
assistance  given  by  Professor  E.  J.  Andrews    is  greatly  appreciated  by 
the  writer. 
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1.  Introduction. 

This  study  is  one  of  a  series  being  conducted  at  the  United  States 
Naval  Postgraduate  School  to  evaluate  the  use  of  automatic  throttle  com- 
pensation systems  in  modern  high  performance  jet  aircraft.     The  need  for 
automatic  throttle  compensation  was  discussed  thoroughly  in  the  previous 
investigations.  [l7  2  J   Briefly,  the  slope  of  the  thrust  required  versus  air- 
speed curve  is  negative  at  the  approach  airspeed  for  modern  jet  aircraft. 
This  means  that  a  loss  of  airspeed  may  rapidly  develop  into  an  uncontrolled 
sink  rate.    Aircraft  with  this  characteristic  are  said  to  have  negative  static 
airspeed  stability.     Positive  airspeed  stability  must  be  provided  by  proper 
use  of  the  throttle.    An  automatic  throttle  compensation  system  relieves 
the  pilot  of  the  task  of  providing  this  airspeed  stability.    With  an  automatic 
throttle  system  keeping  a  constant  airspeed,  the  pilot  can  concentrate  on 
maintaining  his  altitude  or  glide  slope  by  use  of  the  elevator. 

There  are  various  combinations  of  errors  in  airspeed,   pitch  attitude, 
angle  of  attack,  and  normal  acceleration  that  can  be  used  as  inputs  to  an 
automatic  throttle  device.     Direct  use  of  these  feedbacks  and  their  deriva- 
tives and  integrals  provide  the  final  throttle  command  to  the  engine  for  the 
desired  change  in  thrust.    The  system  under  consideration  here  utilizes 
feedback  of  angle  of  attack  and  normal  acceleration.    A  block  diagram  of 
this  system  is  shown  in  Figure  1.    The  transfer  functions,  also  given  in 
Figure  1,  are  those  from  Reference  3. 

2.  Discussion. 

Previous  studies  conducted  at  the  U.  S.  Naval  Postgraduate  School 


considered  flight  in  a  wings  level  attitude.     It  is  the  intent  of  this  study  to 
extend  the  analysis  to  include  turning  flight  and  to  investigate  the  possible 
additional  inputs  to  the  basic  automatic  throttle  system  that  may  be  required. 
The  aircraft  parameters  used  are  for  the  F-8,  an  operational  U.S.  Navy  swept - 
wing  fighter.     Previously,   the  equations  involved  were  simulated  on  an 
analog  computer  and  solutions  obtained  in  this     manner.     For  this  investiga- 
tion,  use  was  made  of  the  1604  Control  Data  Corporation  digital  computer 
at  the  U.  S.   Naval  Postgraduate  School.    A  program  titled  "Runge-Kutta 
Solution  of  Ordinary  Differential  Equations  with  Built-in  Input  and  Output 
Routines"  was  used  to  provide  solutions  for  the  equations. 

Since  the  digital  computer  had  not  been  used  in  References  1  and  2, 
the  response  of  the  basic  airframe  in  a  wings  level  attitude  was  first  exam- 
ined.    This  preliminary  investigation  was  conducted  with  and  without  the 
automatic  throttle  compensator  in  use.    When  the  results  of  this  investiga- 
tion showed  close  agreement  with  other  results,   turning  flight  was  investi- 
gated. 

This  analysis  is  confined  to  the  longitudinal  mode  of  operation  and 
the  linearized  three-degree-of -freedom,   body  axes,   equations  of  motion 
derived  from  small  perturbation  theory  are  utilized.     These  equations  as 
derived  in  Reference  4  are: 

X-Force 

*U  *w  Uo  x  A  T 

^    =    ~&u    +    ~^— gcx    -    g  9    +  ~    AT 

m  m  m 
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Z-Force 


ZJL_«    .    Sv^    +    o    +   ±1 


C*   =    -rr-u  +  -=■  CX    +   e  +  -£-  -7 
mU0  m  mU0      ^ 


Moment 

9=   l~u  +~T_0<  +  "X_0<  +  ire+   b  ^  +  _i~  AT 

The  values  for  the  stability  derivatives  in  the  above  equations  were 
obtained  from  Reference  3  and  are  given  in  Table  I.     These  values  are  for 
the  F-8  aircraft  in  landing  configuration  at  approach  speeds    Several  of  these 
parameters  are  functions  of  Cl  and  Cq  and  would  change  values  if  found 
advisable  to  change  the  aircraft  angle  of  attack  while  in  a  turn.     Utilizing 
the  values  of  the  parameters  and  their  basic  definitions,   they  are  developed 
in  Appendix  I  as  functions  of  Cj^  and  Cjj. 

Curves  of  C^  and  Cjj  versus  angle  of  attack  in  Reference  5  were  used 
to  obtain  the  plot  of  Cj^  versus  Cjj  shown  in  Figure  2.    An  equation  of  this 
plot  is  derived  in  Appendix  II  and  found  to  be 

CD   =    .955    -    2.12(CL)    +    1.44(CL)2 

In  order  to  write  a  program  for  the  digital  computer  a  flow  graph  of 
the  system  to  be  analyzed  is  first  developed.    Such  a  flow  graph  for  the  air- 
frame and  automatic  throttle  combined  is  shown  in  Figure  3»    The  program 
for  the  digital  computer  developed  from  this  flow  graph  is  shown  in  Figure  4. 
Minor  variations  in  this  program  and  the  use  of  initial  conditions  allows  for 
an  analysis  of  the  system  under  a  great  variety  of  conditions. 
3.         Flight  with  Wings  Level. 
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Aircraft  Without  Automatic  Throttle 

By  eliminating  the  change  of  thrust  inputs  to  the  moment  and  X-force 
equations,   the  automatic  throttle  will  be  disconnected  from  the  aircraft 
equations  of  motion.    This  was  easily  effected  by  changing  several  cards 
of  the  program  of  Figure  4.     The  response  of  the  airframe  to  inputs  of  an 
initial  velocity  gust,  an    initial  angle  of  attack,  and  a  step  input  of  elevator 
was  analyzed. 

Figures  5a,    5b,  and  5c  show  the  response  of  the  airframe  velocity, 
angle  of  attack  and  pitch  attitude  respectively  of  an  initial  condition  of 
velocity  equal  to  ten  feet  per  second.     From  this  response,  the  phugoid 
period  is  found  to  be  33  seconds.    This  is  in  close  agreement  with  that 
obtained  in  Reference  1    which  was  33.5  seconds. 

The  airframe  velocity,  angle  of  attack  and  pitch  attitude  response  to 
an  initial  angle  of  attack  of  two  degrees  is  shown  in  Figures  6a,   6b,  and 
6c.     The  short  period  mode  measured  from  Figure  6b  is  6  seconds.    This 
value  also  agrees  with  that  found  in  Reference  1    which  was  6.  06  seconds. 

Figures  7a,   7b,   7c,  and  7d  show  the  airframe  response  to  a  step  input 
of  one  degree  of  elevator.    A  relationship  between  the  angle  of  attack  and 
elevator  can  be  obtained  from  Figure  7b.    A  similar  relationship,   developed 
in  Appendix  III,  agrees  with  that  from  Figure  7b.     This  relationship  is: 
-Y7     =    .46o< 

Aircraft  With  Automatic  Throttle 
The  response  of  the  aircraft  with  the  automatic  throttle  connected  was 
investigated  with  the  disturbance  of  an  initial  velocity  gust  of  ten  feet  per 
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second.     This  disturbance  was  introduced  with  no  elevator  in  the  system  and 
with  the  elevator  function  developed  above  as  a  forcing  function  in  the  system. 
The  aircraft  velocity  and  angle  of  attack  and  the  engine  change  of  thrust  re- 
sponses with  no  elevator  introduced  are  shown  in  Figures  8a,   8b,  and  8c. 
These  same  responses  with  the  elevator  input  is  shown  in  Figures  9a,   9b, 
and  9c.     The  use  of  the  input  of  elevator  as  a  function  of  angle  of  attack 
reduced  the  required  change  of  thrust  by  a  small  amount.     Use  of  this  eleva- 
tor was  made  during  the  rest  of  the  analysis. 

The  velocity,  angle  of  attack  and  change  of  thrust  responses  to  an 
initial  condition  of  two  degrees  to  the  aircraft  pitch  attitude  is  shown  in 
Figures  10a,    10b,  and  10c. 

The  aircraft  response  with  and  without  the  automatic  throttle  as  shown 
in  Figures  5  through  10  agrees  with  the  results  obtained  in  References  1,   2, 
and  3.    This  indicated  that  the  use  of  the  digital  program  simulates  the 
system  to  be  analyzed  as  well  as  the  analog  computer  utilized  in  the  pre- 
vious investigations. 
4.         Turning  Flight  at  Constant  Altitude. 

As  an  aircraft  enters  a  turn  and  maintains  an  angle  of  bank,   more  lift 
must  be  generated  in  order  to  maintain  altitude.    This  is  due  to  the  fact 
that  the  aircraft  lift  is  no  longer  acting  vertically,   but  at  an  angle    <()  ,  while 
the  weight,  of  course,   still  acts  vertically.    The  increased  lift  required  can 
be  generated  by  increasing  the  velocity  while  maintaining  the  same  angle 
of  attack  or  by  increasing  the  angle  of  attack  while  maintaining  the  same 
velocity  or  by  a  combination  of  an  increase  in  velocity  and  angle  of  attack. 
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Figure  1 1  is  a  plot  of  the  thrust  required  for  flight  for  angles  of  bank  up  to 
thirty  degrees.     The  calculations  for  this  plot  are  shown  in  Appendix  IV. 

A  cross -plot  of  the  change  in  thrust  against  angle  of  bank  at  constant 
angle  of  attack  is  plotted  on  Figure  12.    An  equation  of  change  in  thrust  as 
a  function  of  angle  of  bank  was  developed  from  this  curve  and  found  to  be 
AT    =    .05  <|>    +    .855(J>2 

As  an  aircraft  maintains  altitude  in  a  turn  by  generating  more  lift,  the 
normal  acceleration  must  necessarily  increase.    This  increase  in  normal 
acceleration,  n,  as  a  function  of  angle  of  bank  is  easily  found  as  follows: 

cos  <b  =  — — 

T         1  +  n 

1  ■■+  a  «   — br 
COS  $ 

n    =    SEC  <))    -  1 

In  order  to  keep  the  automatic  throttle  from  decreasing  the  thrust 
because  of  this  increase  in  normal  acceleration,  the  base  for  the  normal 
acceleration  loop  must  be  changed.    This  can  be  accomplished  by  putting 
another  input  into  this  loop  or  "biasing"  the  loop.     Figure  13  shows  the 
block  diagram  of  the  system  whith  the  expected  input  changes  as  the  air- 
craft goes  into  a  turn.    The  angle  of  attack  loop  can  also  be  biased  if  it  is 
decided  to  change  the  angle  of  attack  in  a  turn. 

To  investigate  the  use  of  an  automatic  throttle  device  in  turning  flight, 
the  perturbations  introduced  tothejthreedongitudinal  equations  of  motion  due 
to  the  turn  must  be  determined.    An  analysis  of  the  possible  velocity  and 
angle  of  attack  perturbations  is  made  in  Appendix  V.    After  an  aircraft  enters 

14 


a  turn,   if  the  proper  perturbations  are  introduced  and  the  throttle  and  eleva- 
tor are  coordinated  in  such  a  way  as  to  correct  for  these  disturbances,  a 
steady  state  should  be  reached  at  which  the  angle  of  attack,   velocity  and 
thrust  level  out  at  the  new  desired  values. 

It  was  decided  to  make  the  turn  at  a  constant  angle  of  attack,  thereby 
eliminating  any  need  for  biasing  the  angle  of  attack  loop.    This  also  elimi- 
nated the  need  to  make  any  of  the  aircraft  parameters  variable.   Utilizing 
the  system  shown  in  Figure  13,   the  computer  program  was  revised  to  put 
the  aircraft  into  thirty  degrees  angle  of  bank  at  a  rate  of  roll  of  ten;  degrees 
per  second. 

Figures  14a,    14b,  and  14c  show  the  response  of  the  aircraft  to  enter- 
ing this  thirty  degree  angle  of  bank  turn  with  no  perturbations  introduced  to 
the  airframe  due  to  the  turn.    When  the  velocity  perturbation  was  introduced 
alone,  the  response  obtained  was  that  shown  in  Figures  15a,    15b,  and  15c. 
With  both  the  velocity  and  angle  of  attack  perturbations  introduced,  the 
response  was  that  shown  in  Figures  16a,    16b,  and  16c. 

Of  these  three  responses,  the  one  with  only  the  velocity  perturbation 
introduced  produced  the  results  most  closely  agreeing  with  the  desired 
steady  state  as  shown  in  Figure  11.     The  response  with  no  perturbations 
introduced  shows  the  changes  of  velocity,  angle  of  attack  and  thrust  all 
tending  toward  zero  after  90  seconds.    When  both  velocity  and  angle  of 
attack  perturbations  were  used  the  response  of  the  aircraft  shows  the  velocity 
and  angle  of  attack  both  tending  toward  zero  while  the  thrust  has  decreaaedo 
more  than  9000  pounds  and  is  still  decreasing.    The  response  of  the  aircraft 
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with  only  the  velocity  perturbation  introduced  shows  the  velocity  steadying 
at  about  19  feet  per  second,   the  angle  of  attack  approaching  zero  and  the 
thrust  change  steadying  at  about  1450  pounds.     This  is  close  to  the  desired 
steady  state  with  the  undesired  being  that  the  thrust  has  steadied  at  a  value 
about  680  pounds  too  high. 
5.         Conclusions  and  Recommendations. 

During  this  investigation  the  use  of  the  digital  computer  was  found 
very  satisfactory  in  that  originally  putting  the  system  into  a  useable  program 
was  much  easier  than  it  would  have  been  to  wire  an-    analog  board.    Also, 
it  was  not  necessary  to  do  any  scaling  of  the  equations  for  the  digital  com- 
puter, as  it  would  have  been  for  the  analog  computer.     It  was  found  to  be 
quite  easy  to  make  changes  to  the  system  when  the  digital  program  was  used. 

The  difficult  part  of  this  investigation  was  in  simulating  the  aircraft 
making  a  turn.    The  perturbations  introduced  to  the  longitudinal  equations 
of  motion  by  angle  of  bank  are  not  obvious  and  a  completely  satisfactory 
simulation  of  the  turn  was  not  realized  in  this  investigation.    A  more  thorough 
investigation  of  the  effect  of  angle  of  bank  on  the  longitudinal  equations  of 
motion  is  needed  to  reach  a  definite  conclusion  as  to  what  changes  must 
be  made  in  the  basic  automatic  throttle  system  to  handle  a  turn  properly. 

Another  perturbation  that  possibly  should  have  been  included  in  this 
investigation  is  the  rate  of  change  of  pitch  attitude.    Appendix  VI  develops 
a  relationship  between  the  pitch  rate  and  the  normal  acceleration  which  is 
itself  a  function  of  the  angle  of  bank  of  the  aircraft.     Including  this  pertur- 
bation may  simulate  the  aircraft  in  a  turn  better  than  simulations  included 
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in  this  report,     The  perturbations  are  introduced  to  the  airframe  equations  of 

motion  in  the  following  manner. 

X-Force 

*u  /  ,         XwUo  v  XAT 

u    =     ~ K  +  u)    +-— -Jo<t   +  *<)  -g9    +    —     AT 


Z-Force 


Zu ,  %£  '  .        .        Z7 


*  =  mu^(ut  +  u)  +  —  (^t+°<)  +  (et+e)  +  muo  ^ 

Moment 

B-     1li(ut+u)  +  ^-^o<+-^-SK4<>c)  +  -ia-(9t+9) 
M-o  M  aT 

+  T*7  +  ~f  *T 

It  is  recommended  that  any  further  investigations  along  these  lines 
be  continued  as  indicated  in  the  preceding  equations .     It  is  further  recom- 
mended that  the  possibility  of  using  the  hybrid  analog-digital  computer  be 
investigated.    This  system  would  have  the  advantages  of  both  computers  and 
allow  a  great  amount  of  versatility. 


17 


BIBLIOGRAPHY 

1.  Bell,   G.  R.    An  investigation  of  the  Effect  of  Auto-Throttle  Devices 

on  Aircraft  Control  in  the  Carrier  Landing  Approach.     U.  S.  Naval 
Postgraduate  School,    1963  -  B  360. 

2.  Evans,   R.  E.  and  Schuppe,   R.  H.    An  Investigation  of  the  Effect  of 

Simplified  Equations  of  Motion  and  Various  Inputs  on  Automatic- 
Throttle  Devices.     U.  S»  Naval  Postgraduate  School,    1964  -  E  775. 

3.  Anon.    Comparison  of  Two  Automatic -Throttle  Controllers  for  the  F8 

U-2N  Aircraft.     U.   S.  Naval  Air  Development  Center.     Report  No. 
NADC  -  ED  -  LG  292,    1962. 

4.  Etkin,  B.     Dynamics  of  Flight.     John  Wiley  and  Sons,   Inc.,    1959 

5.  Anon.    Criterial  for  Predicting  Landing  Approach  Speed  Based  on  an 

Analog  Computer  Analysis  of  21  Jet  Propelled  Aircraft.     Vought 
Aeronautics.     Report  No.  EOR  -  13202,    1960 


18 


TABLE  I 
F-8  Aircraft  Parameters 
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7      --,  -0.059     l/sec 
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Figure  1 
Block  Diagram  of  Automatic  Throttle 
Airframe  System 
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Figure  3 
Plow  Graoh  of  Automatic  Throttle  Airframe  System 
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. I 

FIGURE  4 

COMPUTER  PROGRAM  OF  AUTOMATIC 
THROTTLE  AIRFRAME  SYSTEM 


PROGRAM  FLY  IT  2 

DIMENSION   X(30)»  XDOT(30)»  C(15) 

C(10)  =  1.0 

PI=4.*ATANF(1.) 

1  CALL  INTEG1     (T»XfXDOT»C) 

THIS  PROGRAM  IS  WITH  AUTO-THROTTLE 
AIRCRAFT  ENTERING  THIRTY  DEGREE  ANGLE 
OF  BANK  TURN 

FOR  LEVEL  FLIGHT  REMOVE  4  NEXT  CARDS 

IF  (T-3.)  2»3»3 

2  PHI=10,*T 
GO  TO  4 

3  PHI=30# 

4  A=COSF(PHI*PI/180. ) 
XU=-.06 

XW=-3.3 
ZU=-. 00113 
ZW=-«426 

ELEV=( .381/.834)*X(2) 

GEE=l./A-l.+(234./32.2)*(X(4)-X<15)  ) 
XDOT ( 1 ) =  XU»X ( 1 ) -32 .2*X ( 3 ) +XW*X ( 2 ) +.00145*X ( 5 ) 
XDOT(2)=ELEV*(-.0  59)+X(4)+ZW*X(2)+ZU»XU) 
X(15)=XDOT(2) 
XDOT(3)=X(4) 

XDOT(4)*ELEV*(-2»5)-.00000455*X<5)-.338*X(4)- 
1.041*X(15)-1»14*X(2)+.000189#X(1) 
XDOT(5)»(l./l«2)*(X(8)+X(6)-X(5) ) 
XDOT(6)=X(7) 

XDOT<7)=-10.*X(6)-11.*X(7)-194000.*GEE 
XDOT(8)*10.*(-X(8)+X(11)+X(10)+X(9) ) 
X(13)=X(2) 

XD0TC9)— X(9)  +  1100«*X(13)*57.3 
XDOT<10)«190.*X(13)*57«3 
XDOT(11)*-2.*X(11)+1140.*X(13>*57.3 
GO  TO  1 
END 
END 


DATA  CARDS  MUST  BE  ADDED  TO  RUN  THIS  PROGRAM 
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Figure  13 

Block  Diagram  of  Revised  Automatic 

Throttle  Airframe  System 
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Aircraft-Automatic  Throttle  Response  to 
Turn  with  Velocity  and 
Angle  of  ^Attack  Perturbations 
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,f  Appendix  I 


If  the  angle  of  attack  is  changed  in  the  turn,  CL 
arid  CD  alscj  change.  The  following  stability  derivatives 
vary  with  Ct  and  CD  and  are  developed  below  as  functions 


of  them,         ■ 


^T   a     -    Cp^oS     r:   -  O.Ofc      *-      C**  .Ml 
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*  -.iscz^.u  +  c*)      yscc 
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Appendix  II 


Figure  2  is  a  plot  of  CL  versus  Cn»     Below  is  a 
derivation  of  a  quadratic  equation  for  Cn  as  a  function 


of  0L. 
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Appendix  III 

Using  the  simplified  relationship  between  elevator 
and  angle  of  attack  given  below  an  equation  for  eleva- 
tor as  a  function  of  angle  can  be  found* 

B 
G^  *  -.381 


^1  /*U*S£  r 


C^7=  -.834 


Cm*, 


(V 


-^  -   .Aloo( 
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Appendix  TV 
Below  are  the  formulae  and  calculations  for  the 
thrust  versus  velocity  at  various  angles  of  attaek. 
The  plot  o^  these  calcutions  is  shown  in  Figure  11. 

U-T«  \f>Ml  SCLaL 

T«  48\o  \bs. 
1.  Keeping  the  same  velocity  in  the  turn. 
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2.  Keeping  the  same  angle  of  attack  In  the  turn. 
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3«     Plying  an  angle  of  attack  corresponding  .to  a  coeffi- 
cient of  lift  of  •67£.   wfcUe  In  the  turn. 
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velocity 
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turn. 
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2.  Keeping  the  same  angle  of  attack  In  the  turn. 

V--^     =   Z4I.4-    fi/sec 

Vtos>  $> 

T=   5\Z.3     U>s. 

3.  Flying  an  angle  of  attack  corresponding  to  a  coeffl< 
elent  of  lift  of   .8821*  while  In  the  turn. 
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C^»  30° 
1.  Keeping  the  same  velocity  In  the  turn. 


»■■'» 


2*  Keeping  the  same  angle  of  attack  in  the  turn. 


3.  Plying  an  angle  of  attack  corresponding  to  a  coeffi« 
cient  of  lift  of  .88214.  while  in  the  turn. 


T^  5UZ5  \ts. 
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Appendix  V 
This   ^ppendix  contains  calculations  of  the  possible 
perturbations  introduced  to  the  velocity  and  the  angle 
of  attack  while. 
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Appendix  VI 

Utilizing  tne  figure  below,  the  equations  following 
ere  easily  derived. 


X-axls 
Y-axls 

Z-axls 
Horizontal  plane 

Turn  rates 


—  COS  (^  =■   Q  S\*\  ^ 
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a  Cv>  ^  0  =   q  cos  <$>  -v  jlL  siv\  ^ 
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Cos     ^  = 


>\*  I 


Combining  the  above  equations  the  following  pitch  rate 
perturbation  can  be  obtained. 
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